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lyser  demonstrated  for  the  first  time. 

•  Optical  visualisation  allows  flow  and 
current  density  to  be  reconciled. 

•  Transition  from  bubbly  to  slug  flow 
associated  with  enhanced  current 
density. 
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A  polymer  electrolyte  membrane  water  electrolyser  (PEMWE)  employing  a  segmented  current  collector 
made  from  a  printed  circuit  board  (PCB)  with  optical  access  to  the  channel  has  been  demonstrated  for  the 
first  time.  The  cell  allows  the  local  current  density,  flow  regime  and  bubble  formation  dynamics  to  be 
studied  in  real  time.  Transition  is  observed  from  a  flow  of  discrete  bubbles  at  the  start  of  the  channel  to 
long  bullet  shaped  bubbles  towards  the  end  of  the  channel  associated  with  a  significant  increase  in  local 
current  density. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Proton  exchange  membrane  water  electrolysers  (PEMWEs)  are 
typically  used  to  produce  hydrogen  from  renewable  energy  sources 
(so  called  ‘green  hydrogen’),  as  they  have  a  millisecond  response 
time  and  are  able  to  withstand  intermittent  electrical  inputs 
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without  degradation.  The  technology  is  therefore  promising  as  a 
component  in  future  national  power  and  transportation  fuel  sys¬ 
tems.  The  PEMWE  technology  is  cited  in  the  European  Hydrogen 
and  Fuel  Cell  programme  as  a  key  process  for  transforming  zero- 
carbon  electricity  sources  into  the  supply  of  zero-carbon 
hydrogen  and  oxygen  for  miscellaneous  end  uses  [1]  and  can  act 
to  accelerate  the  development  of  a  hydrogen  economy  [2]. 

PEMWEs  convert  electrical  energy  to  chemical  energy  by 
applying  a  DC  voltage  across  an  electrochemical  cell.  The  anode  and 
cathode  are  separated  by  a  proton  exchange  membrane  (typically 
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Nafion)  forming  a  zero  gap  electrolyser,  gas  diffusion  layers  (GDL) 
are  used  to  supply  current  and  aid  diffusion,  bipolar  plates  and  end 
plates  hold  the  cell/stack  together.  In  such  a  system,  water  is 
introduced  at  the  anode  where  it  is  oxidised  to  oxygen,  with  the 
production  of  hydrogen  ions  and  electrons,  as  shown  in  Eq.  (1). 
Hydrated  hydrogen  ions  (H+.xH20)  migrate  to  the  cathode  through 
the  membrane  where  they  recombine  electrochemically  with 
electrons  to  produce  hydrogen  (Eq.  (2));  the  overall  reaction  cor¬ 
responds  to  Eq.  (3)  [3-6]. 

Anode:  2H20^02  +  4H+  +  4e~  (1) 

Cathode  :  4H+  +  4e~  ^2H2  (2) 

Overall:  H20->H2  +  1/202  (3) 

PEMWE  technology  has  many  advantages  over  other  electrolysis 
methods  such  as  alkaline  and  solid  oxide  electrolysis  (SOEC)  [7]. 
They  are  capable  of  achieving  current  densities  over  2  A  cm”2, 
reducing  operational  costs  and  potentially  the  overall  cost  of  elec¬ 
trolysis.  The  thin  polymer  electrolyte  membrane  provides  good 
proton  conductivity,  reducing  Ohmic  losses,  and  has  a  low  gas 
crossover  rate,  yielding  hydrogen  with  high  purity  [8].  Further¬ 
more,  PEMWEs  are  able  to  operate  under  high  differential  pres¬ 
sures  (up  to  350  bar)  9  and  have  a  compact  and  simple  design  10]. 

The  first  PEM  electrolyser  was  developed  by  General  Electric 
(USA)  in  1966  for  the  purpose  of  oxygen  generation  [11,12]. 
Currently,  Giner,  Inc.  (USA)  and  Proton  OnSite  (USA)  are  the  two 
major  developers  in  the  USA  and  report  technologies  capable  of 
achieving  system  efficiencies  close  to  70%  (LHV)  [13,14].  However, 
the  cost  of  hydrogen  produced  by  commercially  available  PEMWE 
systems  is  currently  estimated  to  be  $4.8/kg_1-H2  which  exceeds 
the  Department  of  Energy  (DoE)  targets  [15].  ITM  Power  is  the 
largest  developer  of  PEMWEs  in  the  UK  and  a  major  exponent  of  the 
UI<H2Mobility  program,  which  promotes  hydrogen  powered  travel 
in  the  UK.  ITM  Power  produce  both  small-scale  self  pressurising 
PEMWE  systems  (HPac  10  and  HPac  40)  as  well  as  larger  systems  up 
to  a  megawatt  (HFuel).  The  systems  are  able  to  produce  hydrogen 
with  99.99%  purity.  The  HFuel  range  is  fully  compliant  and  self- 
contained  for  refuelling  hydrogen  powered  vehicles  or  direct  in¬ 
jection  into  the  gas  grid  as  a  renewable  energy  storage  system. 

Despite  promising  commercial  developments,  the  challenges  of 
high  cost  (of  which  ~70%  is  materials  cost)  and  the  need  to 
improve  efficiency  persist.  A  large  part  of  the  materials  cost  of  the 
electrolyser  comes  from  the  use  of  noble  metals  as  electrocatalysts, 
as  the  acidity  of  the  solid  electrolyte  precludes  the  use  of  non-noble 
metals.  Typically,  platinum  [16],  carbon  supported  platinum  17], 
carbon  supported  palladium  [18]  and  carbon  supported  platinum / 
palladium  19]  are  used  as  the  cathode  electrocatalyst  with  the 
oxygen  evolution  reaction  (OER)  at  the  anode  typically  employing 
iridium,  ruthenium  (or  their  oxides)  [18,20-22  ,  iridium/ruthe¬ 
nium  oxide,  iridium/tin  oxide  [23]  or  platinum/iridium 
electrocatalysts. 

It  is  important  to  have  a  good  understanding  of  the  distribution 
of  current  in  a  PEMWE  to  optimise  cell  design,  achieve  maximum 
performance  and  longevity.  To  date,  there  are  no  studies  present  in 
the  literature  on  experimental  current  density  distribution  (CDD) 
measurements  of  PEMWE  systems.  However,  a  number  of  methods 
for  measuring  local  current  density  distribution  in  PEM  fuel  cells 
(PEMFC)  have  been  reported. 

Stumper  et  al.  [24]  presents  different  techniques  for  the  deter¬ 
mination  of  CDD  in  operating  fuel  cells:  partial  MEA  approach,  sub¬ 
cell  approach  [25]  and  the  current  distribution  mapping  technique 
that  involves  segmentation  of  the  GDL  and  the  use  of  a  passive 
resistor  network.  Segmented  flow  field  approach  in  PEM  [26,27] 


and  solid  oxide  fuel  cells  (SOFC)  [28  ,  segmented  membrane  elec¬ 
trode  assembly  (MEA)  approach  [29]  and  specially  designed 
segmented  measuring  gasket  approach  [30,31  have  also  been 
investigated  for  CDD  measurements.  Two  other  major  CDD  mea¬ 
surement  types  [32]  are  the  Hall  effect  sensor  method  [33]  and  the 
printed  circuit  board  (PCB)  approach  [34-38]. 

PCB  technology  has  been  shown  to  be  an  effective  way  to  make 
fuel  cells  [39,40]  and  perform  current  mapping,  as  the  approach  has 
the  advantages  of  fast  prototype  cycle  times,  low  cost  manufacture, 
high  design  flexibility  and  utilises  low  weight  composite  materials 
[41].  PCB  technology  was  also  implemented  in  our  previous  work 
on  flow  visualisation  42  . 

This  study  employs  a  segmented  PCB  and  MEA  approach  to 
perform  CDD  measurements  on  a  transparent  single  channel 
PEMWE  such  that  two-phase  flow  features  can  be  linked  to  current 
density  distribution. 

2.  Experimental 

2.1.  Synthesis  of  the  catalyst  coated  membranes 

The  Pt-Ru  ink  was  prepared  using  a  commercial  Pt:Ru  50:50 
atomic  %  black  powder  from  Alfa  Aesar,  UK.  The  ink  was  composed 
of  aqueous  dispersion  of  the  catalyst,  40%  Nafion  solution  and 
18  MQ  cm  deionised  (DI)  water  (Millipore  Milli-Q®  system);  the 
resulting  ink  was  mixed  in  an  ultra-sonic  bath  for  2  h. 

Two  sheets  of  untreated  Toray  carbon  paper,  each  with  an  area 
of  5.1  cm2  (0.6  cm  x  8.5  cm),  were  cut  and  separated  into  8  seg¬ 
ments  using  0.1  cm  strips  of  Kapton®  tape.  The  top  and  bottom 
segments  have  an  area  of  0.72  cm2  (0.6  cm  x  1.20  cm)  and  the  rest 
of  the  segments  have  an  area  of  0.54  cm2  (0.6  cm  x  0.9  cm)  each. 

The  brushing  technique  [8]  was  used  to  deposit  the  catalyst  ink 
on  the  carbon  paper  followed  by  application  of  the  decal  method 
[43]  to  prepare  the  catalyst  coated  membranes  (CCMs).  The  carbon 
sheets  were  painted  with  the  catalyst  ink  after  sonication  and  dried 
at  140  °C  for  30  min.  The  strips  of  Kapton®  tape  were  peeled  off  and 
the  coated  carbon  paper  was  weighed.  The  anode  and  cathode 
loadings  were  targeted  to  be  3.0  mg  cm”2.  The  painted  carbon 
sheets  were  then  sandwiched  on  both  sides  of  a  33  cm2 
(3  cm  x  11  cm)  Nafion®  117  membrane  and  hot-pressed  at  500  psi 
pressure  and  120  °C  for  4  min.  The  carbon  papers  on  either  side  of 
the  MEA  composed  were  removed  and  weighed  again  to  calculate 
the  amount  of  catalyst  transferred  to  the  membrane.  Finally,  the 
CCM  prepared  was  immersed  in  18  MQ  cm  DI  water  to  attain 
complete  hydration.  The  membrane  configuration  and  the  di¬ 
mensions  of  the  segments  are  shown  in  Fig.  1(a). 

2.2.  Design  of  the  single  channel  PEMWE  test  cell 

The  anode  and  cathode  flow  fields  were  machined  from  a  1.6  cm 
thick  PCB  using  a  CNC  machine  (Roland  MDX  40).  Each  plate  con¬ 
sisted  of  a  single  channel  (0.3  cm  x  8.5  cm)  and  the  anode  plate  was 
separated  into  8  current  collection  segments.  Shunt  resistors  with 
2  mfi  resistance  were  connected  to  each  segment  of  the  anode.  The 
copper  layers  (35  pm  thick)  of  both  the  anode  and  cathode  boards 
were  electroplated  with  a  0.25  pm  gold  layer  to  avoid  corrosion  and 
to  reduce  contact  resistance.  The  configuration  of  the  segmented 
anode  PCB  plate  is  shown  in  Fig.  1(b). 

Fig.  2(a)  shows  an  exploded  view  of  the  electrolyser  cell 
describing  the  components  used  and  their  assembly.  A  sheet  of 
untreated  (no  Teflon  content)  Toray  carbon  paper  (TGP_H_120)  was 
used  as  the  cathode  GDL  and  the  anode  side  was  left  unsupported 
(no  GDL)  in  order  to  prevent  lateral  conduction  of  current  between 
individual  segments  of  electrode  and  to  gain  optical  access  to  the 
electrode  surface.  The  MEA  and  the  cathode  GDL  were  positioned 
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between  gasket  material  and  flow  field  plates  (that  also  acted  as  the 
current  collectors)  in  such  a  way  that  one  half  of  the  MEA 
(0.3  cm  x  8.5  cm  strip)  was  in  contact  with  the  plates  and  the  other 
half  was  aligned  within  the  flow  channel.  The  cell  was  held  together 
using  end  plates  made  of  2  mm  thick  Perspex  (with  3  mm  wide 
channels)  that  had  12  bolts  going  through  them  which  were 
tightened  to  a  torque  of  1.2  Nm,  as  shown  in  Fig.  2(b). 

2.3.  Measurement  system  and  procedure 

A  schematic  of  the  testing  arrangement  is  shown  in  Fig.  3.  DI 
water  (18  MQ  cm)  from  the  reservoir  was  pumped  through  the  cell 
using  a  programmable  peristaltic  pump.  The  gasses  produced  were 
sent  through  dehumidification  bottles  and  then  to  vent.  The  voltage 
across  each  of  the  shunt  resistors  on  the  anode  flow-field  plate  was 
measures  using  an  18-bit  data  acquisition  board  (National 


Instruments,  UK)  to  calculate  the  current  generated  from  each 
segment.  The  main  anode  and  cathode  flow-field  plates  were 
connected  to  an  Iviumsts.XRi  high  current  potentiostat,  A11700 
(Alvatek  Ltd,  UK). 

2.4.  Current  density  distribution  (CDD)  measurements 

The  electrolyser  cell  was  polarised  to  a  range  of  operating 
voltages  at  different  water  flow  rates  (1,  3  and  5  ml  min-1).  The 
ratio  of  circulating  water  to  the  amount  of  water  required  to 
service  the  process,  (£)  has  a  value  greater  than  1  (£  >  1) 
throughout  all  voltage-current  conditions  when  operated  at  these 
flow  rates  [42]  to  keep  the  membrane  hydrated  throughout  the 
electrolysis  process.  The  current  densities  for  each  current  col¬ 
lector  were  inferred  by  measuring  voltage  drop  across  each  of  the 
shunt  resistors. 


Fig.  2.  (a)  An  exploded  view  of  the  cell  illustrating  the  cell  components  and  their  assembly,  (b)  Illustration  of  an  assembled  electrolyser  cell  used  for  current  mapping 
measurements. 
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Fig.  3.  Illustration  of  testing  arrangement. 


2.5.  Optical  visualisation  measurements 

Images  of  the  liquid  and  gas  flow  patterns  along  the  channel  and 
gas  bubble  formation  were  viewed  using  a  high  speed  video  system 
(Photron  SA-5)  with  Solarc®  light  source  for  illumination.  Images 
were  obtained  at  a  frame  rate  of  2000  fps  and  at  a  resolution  of 
1024  x  1024.  The  images  were  digitally  evaluated  using  Phantom 
Camera  Control  (PCC,  version  v2.14.727.0.)  software. 

3.  Results 

3.1.  CDD  along  the  electrolyser  channel 

Fig.  4  shows  the  CDD  measurements  across  the  8  segments 
under  ‘Low’  (1.45  V),  ‘Medium’  (1.50  V)  and  ‘High’  ( 1.80  V)  operating 
potentials  at  a  water  flow  rate  of  1  ml  min-1.  It  can  be  seen  that  the 
current  density  distribution  along  the  channel  is  relatively  homo¬ 
geneous  at  low  operating  potentials.  However,  while  increasing  the 
operating  potential  of  the  cell  consistently  leads  to  a  local  current 


Fig.  4.  Current  density  distribution  along  the  flow  channel  for  the  8  segments  at  Low, 
Medium  and  High  operating  potentials  when  operating  at  a  water  flow  rate  of 
1  ml  min-1. 
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Fig.  5.  Current  density  of  the  segments  at  Low,  Medium  and  High  cell  potentials  and 
when  operating  at  a  water  flow  rate  of  (a)  3  ml  min  1  and  (b)  5  ml  min-1. 
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Fig.  6.  Two-phase  flow  profiles  of  the  anode  of  the  electrolyser  cell  when  operating  under  Low,  Medium  and  High  potentials  and  water  flow  rates  of  1  ml  min  \  3  ml  min  1  and 
5  ml  min-1. 


increase  at  each  segment,  there  is  a  significantly  larger  increase  in 
current  towards  the  end  of  the  channel. 

Fig.  5(a)  and  (b)  show  the  current  density  variation  across  the  8 
segments  of  the  cell  when  operating  at  water  flow  rates  of 


Fig.  7.  Close-up  of  two  types  of  bubbles  observed  in  the  electrolyser  flow  channel 
showing:  (a)  diameter  measurement  of  spherical-type  bubble  and  (b)  length  of  slug- 
type  (Taylor)  bubble. 


3  ml  min-1  and  5  ml  min  \  respectively.  The  trend  of  CDD  is  similar 
to  that  observed  in  Fig.  4. 

These  results  demonstrate  that  the  increase  in  water  flow  rate 
decreases  the  current  density  variation  across  the  segments  at  low 
and  medium  operating  potentials.  Furthermore,  it  shows  that  at  a 
fixed  operating  potential,  the  current  densities  across  the  segments 
decrease  with  increasing  water  flow  rate,  this  can  be  explained  by 
the  decrease  in  cell  temperature  at  higher  water  flow  rates,  as 
observed  in  our  previous  study  [42  . 

3.2.  Visual  imaging  of  bubble  size  and  two-phase  flow  patterns 

Studies  have  shown  that  the  flow  regime  in  an  electrolyser  can 
have  a  significant  effect  on  performance  [42  .  Flow  pattern  visual¬ 
isation  and  transition  from  bubbly  to  slug  flow  in  an  electrolyser 
has  been  studied  in  previous  work  by  the  authors  [42].  However, 
there  is  little  understanding  of  how  the  evolution  of  flow  affects  the 
local  current  distribution.  By  performing  combined  current  density 


Table  1 

Mean  diameter  and  length  of  bubbles  (and  their  ratios)  observed  at  Low  Medium 
and  High  operating  potentials  and  water  flow  rates  of  1  ml  min-1,  3  ml  min  1  and 
5  ml  min-1. 


Potential 

1  ml  min  1 

3  ml  min-1 

5  ml  min  1 

Db  (mm) 

Db 

k 

Db  (mm) 

Db 

k 

Db  (mm) 

Db 

lb 

Low 

1.23 

1.02 

1.07 

1.00 

0.87 

1.00 

Medium 

2.57 

1.01 

2.00 

1.01 

1.85 

1.01 

High 

2.93 

0.60 

2.92 

0.65 

2.92 

0.72 
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Fig.  8.  (a)  Cross  section  of  the  electrolyser  channel;  (b)  flow  dynamics  of  a  rising  bubble  in  a  channel  with  lower  void  fraction;  and  (c)  higher  void  fraction. 


distribution  and  optical  visualisation  of  bubble  characteristics  it  is 
possible  to  reconcile  the  two. 

Flow  visualisation  is  capable  of  characterising  the  local  flow 
regime:  bubble  flow  is  characterised  by  uniform  distribution  of 
discrete  bubbles  with  diameter  considerably  smaller  than  the 
channel  width  ( Dt>  <  wch)  in  a  continuous  liquid  phase,  whereas 
slug  flow  is  characterised  by  ‘bullet  shaped’  bubbles  which  have  a 
diameter  almost  equal  to  the  channel  width  (Dt,  ~  wch)  but  with 
elongation  along  the  channel  direction. 

Therefore,  from  first  principles  it  can  be  determined  that  the 
criterion  for  bubbly  flow  is: 


Db_ 

Wch 


<  1 


and  for  slug  flow: 


Wch 


The  spherical  bubbles  in  the  bubbly  region  increase  in  volume  as 
the  flow  transitions  to  the  slug  flow  region,  forming  ‘bullet  shaped’ 
bubbles  with  a  length  k,  giving: 


^<1  (6) 

Lb 

Images  obtained  by  direct  visualisation  of  the  electrolyser  anode 
flow  channel  are  shown  in  Fig.  6.  It  can  be  seen  that  the  regularity  of 
bubble  occurrence  and  their  size  increased  with  increasing  total 


cell  current,  indicating  an  increase  in  oxygen  gas  generation.  It  also 
shows  that  the  bubbles  towards  the  top  part  of  the  cell  (Segments 
6-8)  are  larger  and  more  numerous  (frequent)  than  in  the  bottom 
part  of  the  cell,  this  is  due  to  both  an  accumulation  of  bubbles 
formed  lower  down  in  the  channel  and  a  higher  rate  of  oxygen 
generation  towards  the  top  of  the  channel. 

Furthermore,  the  visualisation  shows  a  decrease  in  bubble  size 
with  increasing  water  flow  rate  at  a  constant  operating  potential 
which  can  be  explained  by  the  natural  decrease  in  void  fraction 
with  increasing  liquid  flow  rate.  To  explore  the  effect  of  increasing 
operating  potential  and  water  flow  rate  on  bubble  size  and  shape, 
representative  images  of  bubbles  under  the  operating  conditions 
considered  were  obtained  at  Segment  7  and  analysed.  It  was 
observed  that  larger  bubbles  had  a  high  tendency  to  coalesce  with 
smaller  bubbles.  In  these  cases,  the  total  volumes  of  the  bubble 
clusters  were  calculated  in  order  to  find  the  equivalent  diameter. 
Close-up  images  of  the  two  generic  types  of  bubbles  seen  in  the 
flow  channel  are  shown  in  Fig.  7.  The  mean  bubble  diameter  and 
length  values  obtained  by  averaging  10  samples  are  given  in  Table  1. 

It  can  be  seen  that  the  ratio  of  bubble  diameter  to  bubble  length 
has  a  value  of  ~  1  at  low  and  medium  potentials  for  all  flow  rates. 
The  mean  bubble  diameter  (at  Segment  7)  increases  from  Low  to 
Medium  operating  potentials  and  approaches  a  critical  value  after 
which  it  deforms  and  starts  growing  vertically,  indicating  a  tran¬ 
sition  from  bubbly  to  slug  flow.  For  Low  and  Medium  conditions, 
increasing  flow  rate  leads  to  a  smaller  average  bubble  diameter  as 
would  be  expected  due  to  the  lower  void  fraction  of  gas  bubbles 
caused  by  higher  liquid  water  flow  rate. 


Sma^er  Taylor  bubble  after 

bubbles  \  coalesce  \ 


Flow 

direction 


Large  bubble 


Fig.  9.  The  sweeping  effect  of  a  Taylor  bubble  rising  up  the  electrolyser  channel. 
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Previous  studies  have  shown  that  transition  from  bubbly  to  slug 
flow  and  bigger  bubbles  (based  on  whole  cell  analysis)  enhance 
mass  transport  (bubble  removal  and  water  access  to  the  electrode) 
in  an  electrolyser  and  hence  results  in  better  performance  [42,44- 
47  .  Fig.  6  shows  that  the  gas  void  fraction  in  the  flow  channels 
increase  with  operating  potential  and  decreasing  water  flow  rate. 
Our  previous  study  discusses  the  influence  of  the  void  fraction  on 
the  mass  transport  coefficient,  which  states  that  it  has  an  inverse 
relationship  [47,48  . 

Based  on  the  correlated  measurements  of  flow  visualisation  and 
current  density  distribution,  a  suggested  explanation  is  proposed  to 
explain  the  link  between  bubble  dynamics  and  localised  current.  A 
rising  bubble  in  an  electrolyser  channel  displaces  water,  creating  a 
transverse  motion  of  water  around  it  (Fig.  8(b)).  Larger  bubbles 
cause  a  greater  displacement  of  water,  which,  at  a  certain  bubble 
diameter,  can  act  to  sweep  bubbles  off  of  the  electrode  surface 
(Fig.  8(c)).  Removal  of  bubbles  that  are  growing  on  the  surface  leads 
to  a  greater  electrode  surface  area  available  for  reaction.  Further¬ 
more,  as  observed  in  Fig.  9,  Taylor  bubbles  (slug  flow)  have  a  ten¬ 
dency  to  coalesce  with  smaller  bubbles  and  remove  them  from  the 
electrode  surface.  Therefore,  it  follows  that  local  transition  from 
bubbly  to  slug  flow  can  be  attributed  to  the  increase  in  mass 
transport  that  results  in  an  increase  in  current  density  along  the 
channel. 

4.  Conclusion 

The  combination  of  localised  current  density  mapping  and  flow 
visualisation  in  a  PEMWE  has  been  demonstrated  for  the  first  time. 
Increased  current  density  towards  the  end  of  the  channel  is  seen  to 
be  correlated  with  a  transition  from  small  discrete  bubbles  to 
bubble  volumes  that  fill  the  channels  (bubbly  to  slug  flow).  The 
results  suggest  that  for  a  given  electrolyser  localised  current  den¬ 
sity  distribution  and  water  feed  flow  rate,  there  will  be  a  charac¬ 
teristic  ‘entrance  length’  along  the  channel  associated  with 
transition  from  bubble  to  slug  flow. 
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